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in ; ABSTRACT 

' An asymmetric surface brightness distribution of the rapidly rotating A7IV-V star, Altair, 

. has been measured by the Navy Prototype Optical Interferometer (NPOI). The observations 

were recorded simultaneously using a triangle of three long baselines of 30m, 37m, and 64m, 
^ ^, on 19 spectral channels, covering the wavelength range of 520nm to 850nm. The outstanding 

characteristics of these observations are (a) high resolution with the minimum fringe spacing 
of l.Tmas, easily resolving the 3-milliarcsecond (mas) stellar disk, and (b) the measurement of 
closure phase which is a sensitive indicator to the asymmetry of the brightness distribution of the 
^3 . source. Uniform disk diameters fit to the measured squared visibility amplitudes confirms the 

Altair's oblate shape due to its rapid rotation. The measured observables of Altair showed two 
features which arc inconsistent with both the uniform-disk and limb-darkened disk models, while 
' the measured observable of the comparison star, Vega, are consistent with the limb-darkened 

disk model. The first feature is that measured squared visibility amplitudes at the first minimum 
do not reach 0.0 but rather remain at «0.02, indicating the existence of a small bright region 
on the stellar disk. The other is that the measured closure phases show non-zero/ 180 degrees 
at all spectral channels, which requires an asymmetric surface brightness distribution. We fitted 
the measured observables to a model with a bright spot on a limb-darkened disk and found the 
observations are well reproduced by a bright spot, which has relative intensity of 4.7%, on a 3.38 
mas limb-darkened stellar disk. Rapid rotation of Altair indicates that this bright region is a 
pole, which is brighter than other part of the star owing to gravity darkening. 

Subject headings: techniques: high angular resolution — techniques: interferometric — stars: individual 
(Altair) — stars: rotation 

1. Introduction pressed as a function of the latitude; the radius 

increases and the brightness decreases from the 
The impact of stellar rotation upon fundamen- p^jg ^o the equator. The first attempt to de- 

tal stellar parameters was originally calculated by termine the stellar oblateness owing to rapid ro- 

von Zeipel (von Zeipel 1924). He showed that the ^^tion was performed with the intensity interfer- 

stellar radius and the surface brightness are ex- 



1 



ometcr at Narrabri (Hanbury Brown ct al. 1967). 
Though they mentioned the possibihty of a rota- 
tionally flattened shape for Altair, their data was 
insufficient to decide the oblateness critically. In 
2001, van Belle et al (2001) observed Altair us- 
ing two baselines of the Palomar Testbed Interfer- 
ometer (PTI, Colavita ct al. 1999). They calcu- 
lated the apparent stellar angular diameters from 
measured squared visibility amplitudes using the 
uniform-disk model and found that the angular 
diameters change with position angle. This was 
the first measurement of stellar oblateness owing 
to rapid rotation. They evaluated the effect of 
surface brightness distribution caused by gravity 
darkening on their angular diameter data. Be- 
cause their simulation showed that the diameters 
are changed only a few %, they neglected the grav- 
ity darkening effect. In 2003, a large oblateness of 
the Be star Achernar was measured with the VLTI 
(Domiciano de Souza et al. 2003). They analyzed 
measured deformation with Roche model which 
takes the surface brightness distribution owing to 
gravity darkening into account, but the deforma- 
tion was too large to be explained within the range 
of the model. Though both the measurement at 
PTI and VLTI demonstrated the determination of 
the stellar parameters such as orientation of the 
rotation axis and the apparent oblateness from the 
squared visibility amplitudes obtained with the in- 
terferometer, the effect of surface brightness dis- 
tribution was not determined. 

Because of the rapid rotation of Altair, where 
apparent rotational velocities range from vsini 
=190km/s (Carpenter et al. 1984) to 250km/s 
(Stoeckley 1968), this star is expected to be grav- 
ity darkened. However, the effects of the stellar 
surface brightness distribution on squared visibil- 
ity amplitudes are small at low-resolution before 
the first minimum. For example, the effect of limb 
darkening on visibilities becomes more significant 
after the first minimum. Consequently, in order 
to discuss the surface brightness distribution of 
rapidly rotating stars, it is indispensable to have 
high-resolution measurements at least over the 
first minimum. In 2001, the NPOI achieved an ar- 
ray configuration (see next section) which for the 
first time allowed for observations with a base- 
line as long as 64 m. With this long baseline and 
visible observing wavelength, the minimum fringe 
spacing becomes as small as 1.7 mas at the sky 



position of Altair. This resolution is high enough 
to observe this «3 mas star well over the first 
minimum. Moreover, the NPOI is equipped with 
a system which enables measurement of closure 
phase in addition to the squared visibility am- 
plitudes since 1996 (Benson et al. 1997). Closure 
phase is a sensitive measure to the asymmetry of 
the brightness distribution of a source. If the grav- 
ity darkened stellar disk is seen neither equator on 
nor pole on, the intensity distribution of the source 
projected on a baseline will be asymmetric and we 
expect that the measurement of the closure phase 
gives us useful information about inclination of 
the star. Consequently, we expect that observa- 
tions of Altair with the NPOI will provide better 
information than earlier observations in order to 
examine the surface brightness distribution of this 
rapidly rotating star. 

Observations arc done for four nights from 25 
May to 1 June in 2001 using three baselines of 30 
m, 37 m, and 64 m of NPOI. We examine the ap- 
parent stellar diameters reduced using a uniform- 
disk model, and discuss the surface brightness dis- 
tribution of this star. 

2. Observations 

The NPOI, located near Flagstaff, AZ, is a long- 
baseline optical interferometer which consists of an 
astrometric subarray (four elements) and an imag- 
ing subarray (six elements) (Armstrong et al. 1998) 
For each baseline, data are collected in thirty-two 
spectral channels cover the observing wavelength 
from 450nm to 850nm. We used three baselines 
consisting of three siderostats; two of which were 
in astrometric subarray, AW, AE, and one of 
which was in imaging subarray, W7. In this paper, 
we state each baseline as follows: OBl(AW-AE), 
OB2(W7-AW), and OB3(AE-W7). The arrange- 
ment and the length of each baseline is shown in 
Tab. 1. 

We observed Altair for four nights from 25-27 

May, and 1 June in 2001 and obtained 19 scans 
(Hummel et al. 1998) of data. Detailed observa- 
tion table is shown in Tab. 2. Interposed with 
the observations of Altair were observations of the 
calibration star, ( Aql, and a comparison star a 
Lyr, Vega. We chose Vega as a comparison star 
because of its similar spectral type, magnitude, 
and angular size to that of Altair. In addition. 
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Table 1 

Arrangement of three baselines used for observation. 



baseline 


station 


E [m] 


N [m] 


Z [m] 


length [m] 


OBI 


AW-AE 


-37.379245 


-2.60124 


0.000314 


37.46964647 


0B2 


W7-AW 


-23.813216 


-17.49076 


-0.080794 


29.54661519 


0B3 


AE-W7 


61.192461 


20.092 


0.08048 


64.40661631 



Table 2: Observation table for Altair. 









OBI 




0B2 




0B3 




date 


scan 


Hour Angle 


P. A. 


UDD 


P. A. 


UDD 


P. A. 


UDD 




no. 


[hr] 


[dcg] 


[mas] 


[dcg] 


[mas] 


[dcg] 


[mas] 


May25 


Pt 81 


-0.470 


182.5 


3.10 


211.6 


3.23 


195.2 


3.16 




Pt 84 


-0.318 


182.8 


3.01 


212.1 


3.16 


195.6 


3.18 




Pt 87 


-0.167 


183.2 


2.99 


212.7 


3.12 


196.0 


3.20 




Pt 90 


-0.017 


183.5 


2.97 


213.4 


3.09 


196.4 


3.19 




Pt 93 


0.132 


183.9 


3.08 


214.1 


3.17 


196.8 


3.15 




Pt 97 


0.629 


185.1 


3.14 


216.7 


3.44 


198.5 


3.18 




Pt 100 


0.827 


185.6 


3.16 


218.0 


3.32 


199.3 


3.18 


May26 


Pt 30 


-0.242 


183.0 


3.04 


212.4 


3.20 


195.8 


3.17 




Pt 35 


0.208 


184.1 


2.91 


214.4 


3.13 


197.1 


3.18 




Pt 38 


0.775 


185.5 


3.20 


217.6 


3.44 


199.1 


3.19 


May27 


Pt 23 


-0.967 


181.4 


3.15 


210.0 


3.25 


194.1 


3.17 




Pt 30 


-0.484 


182.5 


3.14 


211.5 


3.25 


195.2 


3.14 




Pt 35 


-0.146 


183.2 


3.03 


212.8 


3.18 


196.0 


3.21 




Pt 39 


0.142 


183.9 


3.08 


214.1 


3.28 


196.9 


3.18 




Pt 44 


0.475 


184.7 


3.27 


215.8 


3.54 


198.0 


3.17 


JunOl 


Pt 14 


-1.041 


181.2 


3.18 


209.8 


3.26 


194.0 


3.12 




Pt 19 


-0.726 


181.9 


3.09 


210.7 


3.31 


194.6 


3.15 




Pt 24 


-0.410 


182.6 


3.06 


211.8 


3.23 


195.3 


3.17 




Pt 29 


-0.120 


183.3 


3.03 


212.9 


3.23 


196.1 


3.16 



"Position Angle, P. A., is east of north. For OB3, 180 degree are added. UDD is uniform disk diameter. 
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Vega is a good comparison when wc discuss the 
oblateness of Altair because we expect negligi- 
ble deformation from a circularly symmetric disk 
given that its observed rotational velocity is less 
than 20km/s (Freire Ferrero et al. 1983). During 
our observations, the measured squared visibility 
amplitudes of Vega showed dependence on zenith 
angle when the star passed within 10 degrees of 
zenith. No such dependence was found outside this 
region. Consequently, we removed data of Vega 
with zenith angle smaller than 10 degree. Four- 
teen scans of Vega remained after this selection. 
We used all 19 scans of Altair because that star 
doesn't pass the region close to the zenith. We also 
used all data of ( Aql and the number of scans of 
this star during four nights was 24. 

C Aql was used to calibrate measured observ- 
ables of both Altair and Vega. The uniform 
disk diameter of this calibration star is computed 
as O.SOmas from its observed color and magni- 
tude (Mozurkewich et al. 1991). This value is also 
used in other observation (Nordgren et al. 1999). 
Though the size of the star, O.SOmas, is small com- 
pared with Vega and Altair, the star is partly re- 
solved with the longest baseline; measured with 
the full 64.4m baseline, squared visibility ampli- 
tudes become about 0.7 for a perfect interferome- 
ter. Because ( Aql is a rapidly rotating star with 
vsini 320km/s (Royer et al. 2002), we exam- 
ined the possibility that using this star as a circu- 
lar disk with the diameter of O.SOmas biases the 
results of Vega and Altair. We evaluated the ef- 
fect of rotation as 5% changes in diameter and 5% 
asymmetric brightness distribution on the star ow- 
ing to gravity darkening. At the longest baseline, 
we found that these features change the squared 
visibility amplitudes about 6% and phase about 3 
degree at most. The calculated possible errors of 
squared visibility amplitudes and phase were al- 
most the same level with the measurement errors, 
which are described in sec. 2.2. Consequently, 
the effect of the rapid rotation of the calibration 
star is not critical to our result, but we need to be 
careful when we reduce precise physical parame- 
ters directly from measured observables. 

Because the signal to noise ratio of the mea- 
sured visibilities decreases as the wavelength be- 
comes shorter, we used only the 20 reddest of the 
32 spectral channels, thus covering wavelengths 
from 520nm to 850nm. The raw squared visibil- 



ity amplitudes of the calibration star depend on 
channels and baselines. The highest raw was 
0.71 ± 0.04 at the reddest channel of 0B2 {V^ for 
a perfect interferometer on a O.Smas star is 0.96) 
and the lowest raw was 0.13±0.01 at the bluest 
channel of 0B3 {V^ for a perfect interferometer on 
a O.Smas star is 0.6). We also removed the data ob- 
tained at the spectral channel with wavelength, A 
= 633nm, because this channel contains light from 
the NPOI's He-Nc metrology laser. One detector 
for 0B3 was not available during our observations 
and one other one didn't work on 1st June 2001. 
Consequently, the number of data were, 19 spec- 
tral channels x 19 scans=361 squared visibility 
amplitudes on OBI and 0B2, 18 spectral chan- 
nels X 15 scans + 17 spectral channels x 4 scans 
= 338 squared visibility amplitudes on 0B3, triple 
amplitudes, and closure phases. The total number 
of data were 361 x 2 -|- 33S x 3 = 1,738. 

2.1. Detector dead time correction 

Before calibrating the measured observables, we 

corrected the effect of detector dead time. When 
we observe bright stars such as Altair and Vega, 
the effect of detector deadtime td is not negligi- 
ble. The deadtime of the detectors, which are 
used at the NPOl, are approximately 200ns. The 
number of photons counted per unit cycle (500Hz) 
at the detectors of the red spectral channels be- 
come about 800; one photon is detected per ap- 
proximately 2.5 /i sec. In such case, the number 
of photons counted at the detectors becomes less 
than the number of incident photons, and the re- 
sponse of the detector becomes non-linear. This 
is a well known effect and wc calculate the actual 
photon count rate, CRs,, using measured photon 
count rate, CRm, as follows. 



CRa, 



CRt, 



1 — t^CRia 



(1) 



Here, count rate is the number of photons counted 
per one second. 

We calculated the actual visibility using the 
measured visibility as follows, 



IKI = 



IK, 



l-tdCiCx(l-|Kn|2)- 



(2) 



This correction changed the squared visibility am- 
plitudes of Vega up to 20%. In this paper, we treat 
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the detector deadtime correction as described 
above. Triple amplitudes are also corrected based 
on eq. (2). 

2.2. Measurement errors 

The statistical errors of measured quantities 
are calculated and recorded for each scan. How- 
ever, actual measiircmcnt errors arc sometimes 
dominated by long-term errors rather than sta- 
tistical errors within short-term measurement 
(Wittkowski ct al. 2001). Generally, it is not easy 
to evaluate long-term systematic errors. In this 
paper, we calculated the variance of calibrated 
quantities, (Tcai('^j), of the calibration star as long- 
term errors and estimated errors of measured 
squared visibility amplitudes, (Tvs(ii,Aj), as fol- 
lows, 



vs,caU^J 



Eto'lKal(ti,A,)|ViV 



(3) 



Here, ti is the time at the i-th scan, \j is the wave- 
length at the j-th spectral channel, crvs,stat(ii) A^ ) 
is the statistical error, \V{ti,\j)\'^ is the cali- 
brated squared visibility amplitudes of the source, 
|K;ai(ii; Aj )P is the calibrated squared visibility 
amplitudes of the calibration star, and N is the 
number of scans. Wc also estimated measurement 
errors of triple amplitudes as above. For closure 
phases, we used the following equation, 

^%{ti, Aj) = Cf%.st'AtiAj) + C^cp,cal(A)- (4) 

The ratio of errors calculated using eqs. (3), and 
(4) to the statistical errors were 4.8 (for squared 
visibility amplitudes at OBI), 15.4 (0B2), 1.4 
(0B3), 2.7 (for triple amplitudes), 1.1 (closure 
phases). This means that the measurement errors 
of these observations are nearly dominated by the 
long-term errors rather than short-term statistical 
errors. Even if the measurement errors are dom- 
inated by long-term errors, it is confirmed that 
data obtained at NPOI shows no systematic error 
compared with data obtained at other interferom- 
eter (Nordgren,Sudol,and Mozurkewich 2001). 

3. Results 

The squared visibility amplitudes on three base- 
lines, triple amplitudes, and closure phases on 



about 19 channels were obtained simultaneously 
for each scan in our observation. Visibility mea- 
sured at a baseline with projected baseline, Bp, of 
a stellar interferometer, which is the Fourier trans- 
form of the brightness distribution of the source, 
is written as follows. 



J I{x) exp[— ifc-BpxJda; 



(5) 



where I{x) is the brightness distribution projected 
to the baseline, k = 2n/X is the wavenumber. For 
example, visibility of the uniform- disk with the 
angular radius of ruD is written as follows. 



VuD(r-uD) 



2Ji(feBpruD) 
kBpruB 



(6) 



Visibility of the limb-darkened disk with the an- 
gular radius of tld and the linear limb dark- 
ening coefficient of u{k) is written as follows 
(Quirrenbach et al. 1996), 



VLD(r-LD) 



-I- 



-{l-«(fc)} 



2Ji(fcBprLD) 



3 - u{k) ' kBprhD 
6 fiT J3/2{kBpri^r>) 



3 - u{k) 



u{k) 



2 {kBpri^nf/^ ■ 



(7) 



The triple product is the product of the visibil- 
ities on three baselines that form a triangle, 

Vtp = |Vi|exp(-i$i)|F2|exp(-i$2)|1^3|exp(-i$3). 

(8) 

The triple amplitude is the absolute value of the 
triple product, 



\Vtp\ = \Vi\\V2\\V3\, 



(9) 



and the closure phase is the phase of the triple 
product, 

$c = *i + ^2 + $3- (10) 

As we can see from eq. (5), if the brightness distri- 
bution of the source projected to a baseline, I{x), 
is symmetric with x, the imaginary part of the vis- 
ibility is zero and the phase of the visibility takes 
or ±180 degree. If the brightness distribution 
of the source is asymmetric, the imaginary part 
of the visibility becomes non-zero and the phase 
takes non-zero/±180 degree. Generally, it is not 
easy to measure phase of visibility with a ground 
interferometer because of atmospheric turbulence. 
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However, the elosiire phase cancels the effect of at- 
mospheric turbulence and conserves information 
of the source. Consequently, closure phase is a 
useful interfcrometric observable in order to dis- 
cuss asymmetry of brightness distribution of the 
source. 

3.1. Measured squared visibility ampli- 
tudes, triple amplitudes, and closure 
phEises 

Vega was chosen as a comparison star due to its 
similarities in size, spectral type, and magnitude 
with Altair, but yet without the high observed 
rotational velocity which would lead to apparent 
oblateness. Figure 1 shows the squared visibility 
amplitudes and triple products measured during 
the night of 25 May 2001 for Vega. Figure 2 
shows the observables of Altair measured on the 
same night. The dashed lines show the uniform- 
disk model with the fitted angular diameters of 
2ruD = 3.17mas for Altair, and 2ruD = 3.11mas 
for Vega. The solid lines show the limb-darkened 
disk model with the fitted angular diameters of 
2rLD = 3.32mas for Altair and 2rLD = 3.22mas 
for Vega. In order to calculate the limb darkening 
model, we used linear limb darkening coefficients 
calculated by van Hamme (Van Hamme 1993) 
with parameters, Teff=7750K and log g=4.0 for 
Altair, and Tcff=9500K and log g=4.0 for Vega. 

The measured observables of Vega are bet- 
ter fitted with the limb-darkened model (reduced 

= 7.5), than with the uniform-disk model (re- 
duced =17.2). The limb-darkened disk diame- 
ter of Vega measured at PTI was 3.28 ± O.Olmas 
(Ciardi et al. 2001) and the diameters are consis- 
tent with about 2% error. However, Altair was 
not well fitted with either the uniform-disk model 
(reduced x^ = 150) or the limb-darkened model 
(reduced x^ = 154). Altair's large inconsistency 
between the measured observables and models is 
mainly owing to the facts that: (a) the squared 
visibility amplitudes at 0B3 and the triple am- 
plitudes around the first minimum do not de- 
crease to zero, and (b) the closure phases show 
non-zero/180 degree at all spectral channels. In 
addition, the measured squared visibility ampli- 
tudes at OBI are slightly bigger than the circular 
uniform disk model while that at 0B2 is slightly 
smaller than the model. The discrepancy of mea- 
sured squared visibility amplitudes and triple am- 



plitudes from the model around the first minimum 
indicates that there is a small bright component 
on the stellar disk, which is not resolved even with 
the longest baseline. Non-zero/180 degree of clo- 
sure phases means the brightness distribution of 
the source is asymmetric. 

3.2. Apparent Stellar Diameters Reduced 
from the Uniform Disk Model 

Before considering the discrepancy from the 
uniform-disk model, we first examine whether the 
uniform disk angular diameter reduced from the 
squared visibility amplitudes changes with posi- 
tion angle and whether the change is consistent 
with the PTI result. We consider the diameter 
change assuming the elliptical shape of the stellar 
disk. The brightness distribution of an uniform 
ellipse projected to a baseline with position an- 
gle (f) becomes the same with that of the circular 
uniform-disk with angular radius, r{(j)), 

r{(f>) = \Ja? sm^{4> -4>a)+h'^ cos^{4> - (l>o), (11) 

where a is the angular radius at the major axis, b is 
the angular radius at the minor axis, and (po is the 
orientation angle of the ellipse on the sky, where 
(pQ = corresponds to the minor axis pointing to 
the north on the sky. As a result, the squared vis- 
ibility amplitude becomes the same as that of eq. 
(6) with the angular radius, r (</>). Consequently, 
in order to determine the parameters of the ellipse, 
we reduce the data using the following procedure: 
first, calculate the angular diameters at a position 
angle by fitting the squared visibility of each scan 
at each baseline with eq. (6), then compare the 
reduced diameters with eq. (11). Notice r(0) is a 
little different from the intercept of an ellipse at 
an angle (f). 

Though the errors of fitted angular diameters 
are calculated from the errors of visibilities, we 
found that the scatters of the diameters were a few 
times larger than the calculated errors at OBI and 
0B2. Thus we define the errors of the fitted an- 
gular diameters in another way. Considering the 
small change of position angle at each baseline, 4° 
(OBI), 5° (0B2), 8° (0B3) (Tab. 2) during the 
observation, we regard the change in angular di- 
ameter of the ellipse with position angle as small 
compared to the scatter of the fitted angular di- 
ameters. We therefore used the variance of the 
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Fig. 1. — From left to right, squared visibility amplitudes on 0B2, OBI, and 0B3, triple amplitudes, and 
closure phases of Vega measured during 25 May 2001. Six scans were obtained. Dashed lines show the 
uniform-disk model with diameter of S.llmas. Solid lines show the limb-darkening model with diameter 
of 3.22mas. We used linear limb darkening coefficients calculated by van Hamme, with Teff=9500K and 
log g=4.0. The limb- darkening model reproduces the measured squared visibility amplitudes and triple 
amplitudes better than the uniform disk model. 
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Fig. 2. — Prom left to right, squared visibility amplitudes at 0B2, OBI, and 0B3, triple amplitudes and 

closure phases of Altair measured during the night of 25 May 2001. Seven scans were obtained. Dashed lines 
show the uniform-disk model with diameter of 3.17mas. Solid lines show the limb-darkening model with 
diameter of 3.32mas. We used linear limb darkening coefficients calculated by van Hamme, with Teff=7750K 
and log g=4.0. It is significant that the squared visibility amplitudes at 0B3 and the triple amplitudes 
around the first minimum don't decrease to zero, and the closure phases show non-zero/180 degree at all 
spectral channels. 
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fitted angular diameters at each baseline, 0.09mas 
(OBI), O.llmas (0B2), 0.02mas (0B3) as errors of 
the fitted angular diameters. The error of the fit- 
ted angular diameters on 0B3 is smaller compared 
with that of on OBI and on 0B2. We consider this 
fact is owing to that the measured visibilities cover 
the first minimum at 0B3 and the position of the 
first minimum improves the precision of the fitted 
diameters. 

Figure 3 shows fitted angular diameters with 
horizontal axis of position angle. The reduced x'^ 
of 1.5 calculated based on the scatter of the di- 
ameters of Altair, with the circular uniform disk 
model is improved to 0.95 when we allow the 
diameter of Altair to change with position an- 
gle. Fitted parameters with eq. (11) become 
2a = 3.31 ± 0.09mas, 26 = 2.93 ± O.lTmas, = 
35° ± 18°. In our case, the uniform disk model 
estimates the angular diameters of Altair about 
0.14mas smaller than the limb-darkening model. 
Resultant parameters corrected for the limb dark- 
ening effect are consistent with the PTI result 
(van Belle et al. 2001), 2a = 3.461 ± 0.038mas, 
26 = 3.037 ± 0.069mas, ao = 25° ± 9°. Compared 
with Altair, fitted diameters and errors calculated 
from the scatters of the fitted diameters on each 
baseline of Vega are shown in Fig. 3. The diam- 
eter fitted with the circular uniform-disk model 
was, 2ruD = 3.11 ± O.Olmas, x'^ = 1-0 and the 
value of x'^ was not changed when we used the 
elliptical model. We found that Altair is well ex- 
plained by the ellipse rather than the circular disk 
while Vega is well explained as the circular disk. 
The dependence of fitted diameters of Altair on 
the position angle was similar to PTl's. 

3.3. A small bright region on the limb 
darkening disk 

In this section, we examine whether the two fea- 
tures of measured observables of Altair which are 
inconsistent with both the uniform-disk and the 
limb-darkened disk models are reproduced with 
the model of a single bright region on the limb- 
darkened disk. When a bright spot with relative 
intensity of Jp was located at (rp, 0p) in polar 
coordinate on a limb- darkening disk with angular 
radius of rg, visibility is written as follows, 

Kiodei(-^p, r-s, rp, (/>p) = (1 - -fp) yLD(r-s) 

-|-/p exp {-ikBpTp cos(^ - (?ip)} , (12) 



With this model, the visibility amplitude at the 
first minimum becomes /p. We can see the 
value of squared visibility amplitudes around 
the first minimum about 0.02 and expect that 
/p = 0.04 ~ 0.05. Though the existence of bright 
spot changes stellar diameter 2rs, we expect that 
the change of 2rs from 2rLD=3.32mas is small. 
We searched the optimal position of the spot (rp; 
0~ n, Arp=0.01rs, (i>^; 0~360, A^=l [deg]) for 
each set of (/p, 2rs) ranging (/p; 0.034^0.06, 2rs; 
3.25~3.51mas). As a result, we found that the set 
of parameters, (7p=0.047, rs=3.38mas) gives the 
best Xc of 7.3. The x^ value of 7.3 means that 
this model reproduces the measured observables 
at the same level as the limb-darkening model re- 
produced for Vega. Compared with the oblateness 
improved x'^ from 1.5 to 0.95, a spot on the cir- 
cular limb darkened disk improved from 150 
to 7.3. It means that the surface brightness dis- 
tribution is more essential than the oblateness for 
our data. A spot on the circular limb darkening 
disk model gave similar x^ of squared visibility 
amplitudes at OBI and 0B2 to the limb darkened 
elliptical disk model. 

Figure 4 shows the x^ map of (/p, 2rs) and the 
X^ niap of (rp,(/)p) at (/p=0.047, 2rs=3.38mas). 
Though the parameters (/p, 2rs) converge, rp and 
4>-p are not determined independently. These two 
parameters appear in the second term of eq. (12). 
The equation indicates that small change of po- 
sition angle at 0B3 during the observation make 
it difficult to determine rp and independently. 
Consequently, measurement of this star with wider 
range of position angle, or with additional long 
baseline which is located perpendicular to the 
0B3, will be needed to determine the position of 
the bright region accurately. 

4. Discussion 

The measured observables are well repro- 
duced with a bright spot on the circular limb- 
darkened disk model. Because Altair is a single 
star (van Belle et al. 2001) and well known as a 
rapidly rotating star, it is natural to consider this 
spot is a bright pole of the gravity darkened star. 
Van Belle et al. (2001) simulated the effect of the 
gravity darkening as an additional 25% brighter 
region covering 20% of the surface of the star on 
the limb-darkened disk. The relative intensity of 
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Fig. 3. — Uniform-disk diameters of Altair (left) and Vega (right) as a function of baseline position angle. 
Squares show fitted angular diameters at OBI, diamonds show fitted angular diameters at 0B3, and triangles 
show fitted angular diameters at 0B2. Dashed lines show fitted angular diameters of 3.17mas (Altair) and 
3.11mas (Vega). Solid hne shows fitted ellipse with parameters of 2a=3.31mas, 26=2.93mas, and 4>o = 
35° ± 18°. 
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Fig. 4. — Left; map as a function of relative intensity of the bright spot, 7p, and the limb darkened disk 
diameter, 2rs. The minimum value of x^, Xc='''-3i is given for /p = 4.7% and 2rs=3.38mas. Lines show 
2x1, 3Xc! 4Xc- Right; ™ap as a function of position of the bright spot, (rp, ^p), when parameters, 
/s=0.047 and 2rs=3.38mas, are given. 
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Fig. 5. — Left; map as a function of angular diameters at the major and the minor axes of ellipse, 2o, 
26. The minimum value of x^, Xe='''-0' given for 2a=3.77mas and 26=3.29mas. Lines show x^, 2%^, and 
3Xg. Right; map as a function of position of the bright spot, (rp/6, (^p), when parameters 2a=3.77mas 
and 26=3.29mas are given. Lines show 2Xe, 3Xej and 4Xe- 



the additional bright region is 5% which is not 
far off from our result of 4.7%. Considering it as 
a gravity-darkened star, we expect that it is bet- 
ter to take the effect of oblateness in addition to 
the surface brightness distribution into account. 
Then we replaced in eq. (12) with ^01(0), eq. 
(11), and fitted measured observables with the 
model, where a bright spot is on the minor axis of 
a limb-darkened elliptical disk. Fixing the relative 
intensity of the spot, /p =0.047, we searched the 
optimized set of (rp/6, cip) for each set of (2a; 
3.43'-4.25mas, A2a=0.02mas, 2&; 2.96'-3.45mas, 
A26=0.01mas). We found the minimum x^ de- 
creases only a little to 7.0, which is given with the 
set of parameters, 2a=3.77mas, and 26=3.29mas. 
Figure 5 shows the x^ map of (2a, 26) and the x^ 
map of (rp/6, 0p) at (2a=3.77mas, 26=3.29mas). 
It is worth to consider the region where (2a, 26) 
gives smaller x^ than Xc because the number of 
parameters are increased. Two parameters, (r'p/6, 
(/)p), converge better than in the former model 
when parameters (2a, 26) are given. 

Though the decrease of x^ is small, we try 
to consider physical parameters of this star with 
the optimized set of parameters of 2a=3.77mas, 
26=3.29mas, rp/6 = 0.36, 0p = 9°. Assuming a 
rotationally symmetric elliptical stellar body, we 
can calculate the inclination, i = 35° using the 



equation, 

/ 62 ~ 
'^i = \J^--2i^-p- (13) 

In this case, the critical velocity of t;c=430km/s 
(Gray 1976) multiplied by sinz=0.58 becomes 
250km/s. The values wsini=190km/s (Carpen- 
ter et al. 1984), 250km/s (Stoeckley 1968) deter- 
mined from spectroscopy are close to but don't 
break the WcSini. The apparent brightness distri- 
bution with these parameters are shown in Fig. 7. 

5. Conclusion 

We observed Altair with high resolution in- 
cluding the measurement of the triple prod- 
uct using three long baselines of the NPOI for 
four nights. Measured observables indicate the 
asymmetric surface brightness distribution of 
this star; the asymmetry is deduced directly 
from the definition of the visibility and does 
not depend on model. Though the measure- 
ment of the structure of the surface brightness 
distribution has been reported for evolved stars 
(Tuthill,Hanifr,and Baldwin 1997), this is the first 
time that the asymmetric surface brightness dis- 
tribution for a main sequence star has been found 
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Fig. 6. — Measured observables, squared visibility amplitudes, triple amplitudes and close phases of Altair 
with models. Solid lines show a spot on the elliptical limb-darkened disk model. Short dashed lines show a 
spot on the circular limb-darkened disk model. Dashed lines show limb-darkened disk model without spot. 
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from the direct measurement using intcrferomc- 
try. The measured observables are well repro- 
duced with a model of a bright spot of relative 
intensity 4.7% on the limb-darkening disk with 
an angular diameter of 3.38mas. The rapid rota- 
tion of the star indicates that the bright spot is 
a bright pole of gravity-darkened star. Though 
we couldn't determine the position of the pole 
owing to insufficient data, we expect that addi- 
tional observation with sets of long baselines cov- 
ering wider range of position angle will solve this 
problem. More sophisticated gravity and limb 
darkening modeling of a rapidly rotating star 
(Domiciano de Souza et al. 2002) will also help 
the determination of physical parameters when 
we have sufficient multi baselines data. We expect 
these kinds of observations will be realized with 
current and future interferometers and that study 
on rapidly rotating stars will progress with further 
observation by interferometers. 
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porting the use of data reduction system, OYS- 
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ter is a joint project of the Naval Research Lab 
and the US Naval Observatory in cooperation with 
Lowell Observatory, and is funded by the Office 
of Naval Research and the Oceanographer of the 
Navy. 
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Fig. 7. — A possible solution of a bright spot on 
the elliptical limb darkening disk model of Altair 
projected onto the sky with parameters; angular 
diameter at the major axis 2a = 3.77mas, angular 
diameter at the minor axis 2b = 3.29mas, orienta- 
tion of the rotational axis projected onto the sky 
from the north pole 4>o =9°, rcaltivc intensity of 
the bright spot Ip = 4.7%, inclination i = 35°. 



This 2-column preprint was prepared with the AAS lAlJ^X 
macros v5.2. 
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